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ABSTRACT 

We model the evolution of galaxy clustering through cosmic time to investigate the nature of the 
power- law shape of £(r), the galaxy two-point correlation function. While £(r) on large scales is set 
by primordial fluctuations, departures from a power law are governed by galaxy pair counts on small 
scales, subject to non-linear dynamics. We assume that galaxies reside within dark matter halos and 
subhalos. Therefore, the shape of the correlation function on small scales depends on the amount of 
halo substructure. We use a semi-analytic substructure evolution model to study subhalo populations 
within host halos. We find that tidal mass loss and, to a lesser extent, dynamical friction dramatically 
deplete the number of subhalos within larger host halos over time, resulting in a ~ 90% reduction by 
z = compared to the number of distinct mergers that occur during the assembly of a host halo. We 
show that these non-linear processes resulting in this depletion are essential for achieving a power-law 
£(r). We investigate how the shape of £(r) depends on subhalo mass (or luminosity) and redshift. 
We find that £(r) breaks from a power law at high masses, implying that only galaxies of luminosities 
< L* should exhibit power-law clustering. Moreover, we demonstrate that £(r) evolves from being 
far from a power law at high redshift, toward a near power-law shape at z = 0. We argue that 
tf(r) will once again evolve away from a power law in the future. This is in large part caused by the 
evolving competition between the accretion and destruction rates of subhalos over time, which happen 
to strike just the right balance at z « 0. We then investigate the conditions required for £(r) to be 
a power law in a general context. We use the halo model along with simple parametrizations of the 
halo occupation distribution (HOD) to probe galaxy occupation at various masses and redshifts. We 
show that key ingredients determining the shape of £(r) are the fraction of galaxies that are satellites, 
the relative difference in mass between the halos of isolated galaxies and halos that contain a single 
satellite on average, and the rareness of halos that host galaxies. These pieces are intertwined and we 
find no simple, universal rule for which a power-law £(r) will occur. However, we do show that the 
physics responsible for setting the galaxy content of halos do not care about the conditions needed 
to achieve a power law £(r) and these conditions are met only in a narrow mass and redshift range. 
We conclude that the power-law nature of £(r) for L* and fainter galaxy samples at low redshift is a 
cosmic coincidence. 

Subject headings: cosmology: theory — dark matter — galaxies: halos — galaxies: structure - 
large-scale structure of universe 
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1. INTRODUCTION 

The two-point correlation function of galaxies was 
measured four decades ago and found to be consistent 
with a £(r) oc r~ 2 power law (iTotsuii fc Kiharal 119691: 



Peebles 1973; Haus er fc Peebleslll973t iPeebles fe Hauserl 
1974 IPeebles! Il974| ). Since that time, successively 
large r galaxy redshif t surv ey s (e.g.. |Huchra et all 
19831: IdaCosta et all 119881: SantiagcT gLaLl 119951: 



Shectman et al . 1996; S aunders et al.ll2000l : iColless et all 
20011: [York et al.ll2000t ) have mapped the distribution of 
galaxies with ever increasing precision and confirmed 
correlation functions consiste nt with power laws over 



Marzke et all 11995 


; iHermit et all 19961 iTucker et all 


19971: Dine: et all 


1998. 120021: iNorberg et all 12002: 


Zehavi et al. 20021 The scales on which a single power- 



law description is valid span a range from large regions 
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exhibiting mild density fluctuations (r > 10 Mpc), to 
smaller regions with large density fluctuations experi- 
encing rapid non-linear evolution (r ~ 1 — 10 Mpc), 
to collapsed and virialized galaxy groups and clusters 
(r < 1 Mpc). It has long been noted that the lack 
of any feature delineatin g the transitions among these 
scales is surprising (e.g.. IPeebles! Il974t iGott fe Turne: 
1979t lHamilton fc Tegmarkl 1200% IMasiedi et all 1200 



Li fc White! 120091 ). This is especially true given that the 



matter correlation function in the now well-established 
concordance cosmological model differs significantly 
from a power law. In this paper, we return to this long- 
standing problem and address the origin of a power-law 
galaxy correlation function in the context of our modern 
paradigm for the growth of cosmic structure. 

This conundrum can be refined within the con- 
temporary framework in which galaxies live within 
virialized halos of dar k matter (jWhite fc Reesl Il978t 
iBlumenthal et al. 1984). In such a model, galaxy clus- 
tering statistics can be modeled as a combination of dark 
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matter halo properties and a halo occupation distribution 
(HOD) that specifies how galaxies occupy their host ha- 
los fe.g.JPeacock fc SimthlboOOtlScoccimarro et al.ll2001t 
IBerlind fc Weinberd 12003: ICoorav fc Shethl 120021 . In 
this halo model approach, the galaxy correlation func- 
tion is a sum of two terms: On small scales, pairs of 
galaxies reside in the same host dark matter halo (the 
"one-halo" term), whereas on large scales, the individual 
galaxies of a pair reside in distinct halos (the "two-halo" 
term). These two terms depend on the HOD in differ- 
ent ways, requiring delicate tuning in order to spaw n 
an unbroken power law (e.g.. IBerlind fc Weinber j l2002). 
Consequently, a feature in £ (r) at scales corresponding to 
the radii of the typical, viralized halos that host luminous 
galaxies is expected. 

In a dramatic success for the halo model, IZehavi et al.l 
(2004) first detected a statistically-significant departure 
from a power law due to the high precision measurements 
of the Sloan Digital Sky Survey, and demonstrated that 
the halo model pro vides an acceptable fit to the data. 
IZehavi et al.l ([2005) confirmed this result, adding that 
power-law dep artures grow stro nger with galaxy lumi- 
nosity (see also lBlake etaT1l2008l : iRoss et al.ll2010f ) . £(r) 
has since been shown to devi a te from a power l aw at high 
redsh ifts (lOuchi et al.l |2005t iLee et alj[2006t ICoil et al.l 
2006; IWake et al.ll2011[ L Nevertheless, it remains a fact 
that deviations from a power law at low redshifts are 
small and the galaxy correlation function is roughly a 
power law over an enormous range of galaxy-galaxy sep- 
arations. Deviations have been revealed only through 
ambitious observational efforts. 

Halos are known to be repl ete with self-bound 
structures, dubbed "subhalos" (jGhigna et al.l 119981 : 
iKlvpin et ail Il999t iMoore et all [1999J), and both halos 
and subhalos are thought to be the natural sites of galaxy 
formation. Subhalos were isolated halos in their own 
right, hosting distinct galaxies before merging into a 
larger group or cluster halo 1 . Remarkably, the clustering 
of host halos along with their associa ted subhalos is very 
similar to that of observed galaxies ([K ravtso v fc Klvpinl 
119991: ICohn et al.lll999l : iKravtsov et al.ll2004af ). suggest- 
ing a simple correspondence of galaxies with host ha- 
los and subhalos. This was clearly demonstrated by 
IConrov et al.l (|2006| ) who compared the correlation func- 
tions of hosts and subhalos to that of galaxies over a 
broad range of luminosities and redshifts (z ~ — 4), 
finding excellent agreement. These results indicate that 
an understanding of the physics governing the subhalo 
populations within host halos may provide insight into 
the physics of galaxy clustering and the near power-law 
form of the galaxy two-point correlation function. 

In this paper, we examine the causes of the observed 
power-law correlation function by studying the mergers, 
survival, and/or destruction of dark matter subhalos. 
Our focus in this paper is on the gross features of the 
galaxy two-point function and not on detailed compar- 
isons to specific data sets. We explore more sophisticated 
galaxy-halo assignments and statistical comparisons with 
data in a forthcoming follow-up study (Watson et al. in 

1 Satellites or subhalos are used throughout the paper to refer 
to self-bound entities lying within the virial radius of a larger halo. 
Those that do not lie within a larger system are designated as 
centrals, host halos or simply hosts. 



prep.). 

We argue that the nearly power-law, low-redshift 
galaxy correlation function is a coincidence. The corre- 
lation function of common L < L» galaxies evolves from 
relatively strong small-scale clustering at early times, 
through a power-law at the present epoch, and most 
likely toward relatively weak small-scale clustering in the 
future. The origin of the present-day power law, in turn, 
relies on the tuning of several disconnected ingredients, 
at least three of which are: the normalization of primor- 
dial density fluctuations determined by early Universe 
physics; a halo mass scale for efficient galaxy formation 
determined largely by atomic physics, stellar physics, and 
the physics of compact objects; and relative abundances 
of baryonic matter, dark matter, and dark energy in the 
Universe. 

Our paper is organized as follows. In § 12.11 we review 
the halo model and restate the problem in terms of this 
framework. In § [3] we give an overview of our primary 
modeling technique. In §[4] we investigate the individual 
roles of merging, dynamical friction, and mass loss in 
shaping the halo occupation statistics of subhalos, as well 
as the resulting halo correlation function. In §[5] we show 
how £(r) depends on host halo mass and redshift. In §[B] 
we explore a standard parametrization of the HOD to see 
what is required to get a power-law £(r), and we predict 
the masses and redshifts at which a power-law £(r) can 
be constructed. In §[7] we give a summary of our results 
and our primary conclusions. Throughout this paper, 
we work within the standard, vacuum-dominated, cold 
dark matter (ACDM) cosmological model with Q, m = 0.3, 
n A = 0.7, Q h = 0.04, h = 0.7, <j 8 = 0.9, and n s = 
1.0. These values differ slightly form the WMAP best- 
fit values, however this has little effect on our general 
results and was chosen in order to compare to previous 
work that used similar cosmological models. 

2. A MODERN RESTATEMENT OF THE PROBLEM IN 
HALO MODEL LANGUAGE 

Though the observed galaxy correlation function is 
nearly a power law, the matter correlation function pre- 
dicted by the concordance cosmological model is not. 
This is evident in Figure [TJ On scales corresponding to 
collapsed objects, the dark matter correlation function 
exceeds the values that would be obtained by extrapolat- 
ing the larger-scale power law to small scales. However, 
galaxies are biased with respect to dark matter in such 
a way as to counteract this excess. We can examine this 
discrepancy in terms of the halo model. If the reader is 
familiar with the halo model formalism, he or she may 
wish to skip to § 12.21 

2.1. Halo Model Basics 

Assuming that all galaxies live within virialized dark 
matter halos, the galaxies comprising any pair can come 
either from within the same halo (the one-halo term) 
or from two separate halos (the two-halo term). The 
correlation function is then given as the sum of these 
two terms 

£(r)=£(r) lh +£(r) 2h + l, (1) 

(e.g., ICoorav fc She th 2002; for this particular form of 
the equation see [Zheng 2004). The probability distribu- 
tion P{N\M) that a halo of mass M contains N galaxies 
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Fig. 1. — Correlation function of galaxies compared to dark 
matter. Points show the correlation function of galaxies from the 
APM sur vey, estimated from deprojecting the angular correlation 
function (Maddox ct al. 1990: Baugh 1996). The curve shows the 
correlation function of dark matter measur ed from the LCDM G IF 
simulation run by the Virgo collaboration (Jenkins ct al. 1998). 



together with the spatial distribution of galaxies within 
their host halos constitute the halo occupation distribu- 
tion (HOD). We denote the first and second moments of 
P(N\M) at a specific mass M as (N) M and (N(N-1)) M , 
respectively. The one-halo term can be computed by 
counting the average number of galaxy pairs of a given 
separation in a common halo and a veraging over all ha- 
los. W e write the one-halo term as ([Berlind fc Weinberg) 
2002) 



l + £(r) 



lh 



(N(N — 



1 S ))m 



dM^p- (2) 
dM y ' 

F(r\M) 



where dn/dM is the halo mass function, (N(N — 1))m/2 
is the mean number of galaxy pairs within a halo of mass 
M, and F{r\M) is the distribution of separations be- 
tween these pairs 2 . If the average spatial distribution 
of galaxies within their host halos is X(r\M), then the 
pair separation distribution F(r\M) is the convolution 
of X(r\M ) with itself. The quantity n g is the mean den- 
sity of galaxies in the Universe, 



M ■ 



(3) 



Motivated by theoretical considerations (e.g., 
Berlind et afl l2003t iKravtsov etafl l2004at IZheng et all 



2005), the HOD of galaxies is usually considered sepa- 



rately for central galaxies that live near the centers of 
their host halos and satellite galaxies that orbit within 
the host halo potential. Each halo above some mass 

2 This notation is s lightly different from that used in 
Berlind & Weinberg (2002), in which F(r) denoted the cumulative 
pair distribution. 



threshold should contain one central galaxy and possibly 
one or more satellites, depending on the host mass 
and the HOD. In this framework it is useful to con- 
sider contributions to the one-halo term separately for 
central-satellite and satellite-satelli te pairs. Therefore, 
we r ewrite the one-halo term as (jBerlind fc Weinbergl 
[200l 



in . 



1 



2nr 2 nl 



(N cen N sat )MF cs (r\M)- 



dM 



(iV sat (jV sa t - 1)) M 



(4) 



F ss {r\M) 



where (N cen N aat ) M and (iV sat (7V sat - l)) M /2 are the 
mean number of central-satellite and satellite-satellite 
pairs in hosts of mass M, and F cs (r\M) and F ss (r\M) 
are the pair separation distributions of central-satellite 
and satellite-satellite pairs, respectively. If the central 
galaxies always reside very close to the center of the host 
halo and the average distribution of satellite positions 
within the host halo is A s (r|M), then F cs — X s (r\M) and 
F ss (r\M) is the convolution of X s (r\M) with itself. In 
practical cases there is at most one central galaxy and 
satellites are only present in halos with a central, so that 
(N cen N sa t) M = (N sat )M. The total fraction of galaxies 
that are satellites in a sample is then 



/sat Tig 



JdM^(N sat ) M 



J dM^((N cen ) M + (N sat ) M y 



(5) 



The satellite fraction, / sat , will prove an important quan- 
tity in determining the shape of the galaxy correlation 
function. 

On scales significantly larger than individual halos, the 
two-halo term dominates the clusterin g strength. It is 
most simply written in Fourier space as (jCoorav fc Shethl 
2002): for this p articular form of the equation see 
Tinker et al.l[2005h 



P 2h (k) = P m (fe) 



(6) 



x b h (M,r)~X(k\M) 



where P m (fc) is the matter power spectrum, bh{M,r) 
is a (possibly scale-dependent) halo bias function, and 
X{k\M) is the Fourier transform of the spatial number 
density of galaxies within their host halos. We can invert 
the Fourier transform of the two-halo power spectrum to 
recover the two-halo term of the correlation function. In 
the limit that the galaxy pair separation is larger than 
any halo of interest, the two- halo term becomes 



e\r), 



: ' bh (M,r)(N) M ^dM 



-1 2 



£m, (7) 



■K 6. 



where £ m (r) is the matter correlation function. Equa- 
tion ([7]) explicitly shows that the large-scale galaxy cor- 
relation function is essentially the halo correlation func- 
tion except halos of different masses are weighted by 
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(N)m- The galaxy bias describing the relative clustering 
of galaxies to dark matter b g = \J £ /£ m is the quantity 
in square brackets in Equation ([7]). 

2.2. The Battle of the 1- Halo and 2- Halo Terms 

iBerlind fc Weinberg! (j2002j ) showed that maintaining a 
power-law correlation function requires a careful balance 
between the one-halo and two-halo terms and is thus 
quite difficult to achieve. This is because the one-halo 
term generally changes by a larger amount than the two- 
halo term in response to changes to the HOD. A close 
examination of Equations ([3]), ((U), ([6]) and (J7J) reveals 
why this is the case. 

Consider first the two-halo term as it is the simplest. 
On large scales, the two-halo term is just a weighted av- 
erage of the clustering of host halos. For simplicity, as- 
sume (albeit incorrectly) the halo bias to be a constant 
function of halo mass. Increasing (N)m increases both 
the number of two-halo pairs at a given separation (the 
square of the integral in Eqs. [S] and [7]) and the num- 
ber of random pairs rig/2, by the same amount. The 
reason the two-halo term is at all sensitive to the HOD 
is that the bias of halos does depend on mass and so 
changing the relative number of galaxies in high-mass 
vs. low-mass halos changes the weight in the average of 
the halo bias in Equation ([7]). For example, assigning 
a large number of satellite galaxies to high-mass halos 
increases £ 2h (r) by weighting highly-biased, high-mass 
halos more heavily. The possible range in the ampli- 
tude of the two-halo term is limited by the variation of 
the halo bias function 6^(M), within the mass range 
relevant to galaxies, 10 11 < M/M Q < 10 15 . At low 
masses, the halo bias is bh ~ 0.65 while, in the cluster 
regime (M ~ 10 14 fe ~ 1 M w ), it grows to values of bh ~ 2 
(Tink er ;t al.l l2005h . Bias continues to grow with mass, 
but more massive halos are rare and do not contribute 
much to the weighted average because dn/dM is minis- 
cule. The two-halo term scales like the square of the 
average bias b g in Equation ((7} , so the possible dynamic 
range £ 2h can display is, at most, a factor of ~ 9 and 
is usually significantly smaller. Simply put, the two-halo 
term depends weakly on the HOD because on large scales 
it is not possible to make galaxies significantly more or 
less clustered than the host halos they occupy. 

On small scales, the one-halo term dominates and the 
situation is different. The number of galaxy pairs within 
an individual halo scales with (N(N — 1))m while the 
number of random pairs scales with n 2 ., or (N) 2 M . It is 
instructive to break the HOD into central and satellite 
galaxies. In the regime where there is one central galaxy 
per halo, the mean number of central-satellite pairs is 
(N cen N sat ) M = (iV Ba t)Af, whereas the mean number of 
satellite-satellite pairs is (N sat (N sat — 1))m/2. Assum- 
ing a Po isson distribution for the number of satellite 
galaxies (jKravtsov et al.l l2004aD . (N sa . t (N sa . t - 1)) M = 
(N sa ,t)\f The mean number of random pairs scales like 
(1 + (iVsat))^. In the limit (N sat )M > 1, the number of 
satellite-satellite pairs dominates the number of central- 
satellite pairs, but in this limit both the number of one- 
halo pairs and the square of the mean galaxy number 
density scale as (N sa t}%[ so the one-halo term saturates 
to a maximum value and is insensitive to the number of 
satellite galaxies per halo. 



In most practical cases, the fraction of satellite galax- 
ies in an observational sample is / sa t < 0.25, so samples 
tend to be dominated by halos with satellite galaxy pop- 
ulations in the opposite limit, (N sa ,t)M 1. This is 
due to the fact that very massive host halos are rare, so 
halos with (N sat ) M > 1 are rare. With (N sat )M < 1, 
the central-satellite term dominates and the number of 
such pairs scales as (N sat )M while the mean number den- 
sity n g is approximately constant. Examination of Equa- 
tions (HJ and <j5j) reveals that in this regime £ lh scales in 
proportion to the fraction of satellite galaxies and in in- 
verse proportion to the number of host halos. Host halo 
mass is largely fixed by requiring the galaxies in any sam- 
ple to have an appropriate average number density (this 
is why rare galaxies exhibit strong small-scale cluster- 
ing). Therefore, the one-halo term describing any given 
sample varies approximately linearly with (N sa .t)M until 
(iVgat) > 1, at which point it saturates. It is interesting 
that nearly all the sensitivity of the correlation function 
to the HOD comes from central-satellite galaxy pairs in 
host halos where satellite galaxies are uncommon! 

In this work, we aim to understand the origin of the 
nearly power-law galaxy correlation function. The rele- 
vant question is why is it that the number of galaxies (or 
satellite galaxies to be more specific) per halo is set just 
so that the one-halo and two-halo terms in the galaxy 
correlation function match smoothly, leaving only small 
deviations from a single power law over several orders of 
magnitude in scale? We confront this problem by study- 
ing the properties and evolution of subhalo populations. 
We now turn to some of the details of our modeling meth- 
ods. 

3. OVERVIEW OF HALO SUBSTRUCTURE MODELING 

Our approach is to study the evolution of subha- 
los within virialized host halos as a method to under- 
stand satellite galaxies and, in turn, the evolution of 
galaxy clustering. We focus our attention on the relative 
strengths of small-scale and large-scale clustering. We 
study subhalo p opulations using the approximate semi- 
analytic model of lZentner et al.l (|2005l hereafter Z05). In 
this section, we briefly review the fundamental aspects of 
the model that are of immediate relevance and we refer 
the reader to Z05 for details and validation. The subhalo 
model is based on lZentner fc Bullock! (12003D and i s simi- 
lar to the in dependent models of | Tavlor fc Babull (2004, 
l2005alTbh and lPenarrubia fc Benson! (|2005f ) . while sharing 
many features with other approximate treatments of halo 
substructure (lOguri fc Leel 120041: Ivan den Bosch et all 
20051 iFaltenbacher fc Mathewsll2005l; iPurcell et alj|2007t 
Giocoli et al 1120081 \200$ T~ 

Semi-analytic models are an approximation to the cal- 
culations of large A-body simulations, yet such mod- 
els offer many advantages: (1) semi-analytic calculations 
are computationally inexpensive; (2) they have no in- 
herent resolution limits; (3) they enable the statistical 
study of subhalos within very large numbers of host ha- 
los; (4) they allow the growth and mass-loss histories 
of particular subhalos to be tracked without significant 
post-processing and analysis; (5) they make studies of 
model parameter space tractable; and (6) semi-analytic 
models facilitate parsing complex physical phenomena 
so that the relative importance of different physical ef- 
fects may be understood. Our goal is to quantify the 
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relative importance of merging, which increases subhalo 
abundances, and dynamical friction and mass loss, which 
decrease subhalo abundances. We also aim to explore 
predictions for subhalo populations and galaxy correla- 
tion functions from high redshift to several Hubble times 
in the future. Z05 extensively tested the model we use in 
this paper and showed that the model produces subhalo 
mass functions, occupation statistics, and radial distri- 
butions within hosts that are in good agreement with a 
number of high-resol ution iV-body simulation s (see the 
recent comparison in iKoushiappas et al.ll2010l as well). 

The analytic model proceeds in several steps. For a 
host halo of a given mass M, observed at a given red- 
shift z, we generate a halo merger tree using the mass- 
cons erving implementation of the excursion set formal- 
ism dBond et al.lll99lULacey fc Cold fl993ll 19941) dev el- 
oped by ISomerville fe Kolattl (|1999l . see lZentnerll2007l for 
a review) . This yields a complete history of the masses 
and redshifts of all halos that merged to form the fi- 
nal, target halo of mass M at redshift z. The host 
halo is the largest halo at each point in the merger 
tre e. We model the de nsity distributions of all halos 
as iNavarro et al.l (|1997l hereafter NFW) profiles with 
concentrat ions determined by th eir merger histories ac- 
cording to IWechsler et al.l (|2002D . At the time of each 
merger, we assign the subhalo initial orbital parame- 
ters drawn from d istributions m easured in N-body sim- 
ulations (Z05, see iBensonl 120051 for similar formalisms). 
We then integrate each subhalo orbit within the host 
halo gravitational field, taking into account dynamical 
friction and mass loss. We e stimate dynami c al fric tion 
with an upda ted form of the IChandrasekharl (l!943[) ap- 
proxi mation (|Hashimoto et all 120031 : iZentner fc Bullock! 
2003j), account for internal heating so that scaling re- 
lations des cribing the internal structures of subhalos 
are obeyed (jHayashi et aD2003l : iKazantzidis et aLll2004t 
iKravtsov et al.l[2004bl ). and allow for loss of material be- 



yond the tidal radius on a timescale comparable to the 
local dynamical time. The details of each ingredient are 
given in Z05. 

The correlation function of halos and subhalos and 
their associated galaxies is sensitive to the abundance 
of subhalos that survive both possible mergers with the 
central, host galaxy due to dynamical friction as well 
as mass loss and thus remain as distinct objects in or- 
bit within their host halos with their galaxies intact. 
Therefore, it is necessary to specify conditions under 
which the galaxy within a subhalo may be "destroyed" 
and removed from our samples. In this work, we con- 
sider the clustering of mass-threshold samples of halos 
and subhalos as a proxy for luminosity-threshold sam- 
ples of galaxies, so significant mass loss will lead to a 
galaxy that is either destroyed or dropped out of our 
sample. We assume such a scaling between halo mass 
and galaxy luminosity solely for the sake of simplicity. 
Our primary points are qualitative in nature, but we 
note that this is similar to ot her schemes that have de- 
scribed data successfully (e.g.. IKravtsov fc Klypinl 1 19991: 
Colin et al.lll999l: IKravtsov et al.ll2004aHTasitsiomi et al.l 



20041: IConrov et al.ll200a) and our calculations with sim- 
ilar, but more sophisticated assignments do not alter any 
of our basic results or conclusions. In rare cases, subhalos 
may survive close encounters with the center of their host 
halo potentials. We remove all subhalos that have orbital 



apocenters r apo < 5 kpc. This choice is physically mo- 
tivated because the galaxies within such subhalos would 
likely have merged with the central galaxy, or at least be 
observationally indistinguishable from the central galaxy. 
This choice is relatively conservative in that galaxies on 
larger orbits would also likely be influenced and it only 
affects the results of calculations in which tidal mass loss 
is not permitted (sec below). The net result of evolving 
orbits for each subhalo in the merger tree is a catalog of 
all surviving subhalos in the final host halo at the time 
of observation. In some cases, a halo that merges into 
a larger host contains subhalos of its own. These subs- 
of-subs are only abundant inside very large host masses 
and are present in our model. 

One of our aims is to study the individual roles of halo 
merging, dynamical friction, and mass loss on the cluster- 
ing of halos. Therefore, we compute subhalo populations 
in four different sets of circumstances: 

No Effects - a "bare-bones" model that does not al- 
low satellite galaxies to be modified by dynamical 
friction or mass loss. In this case, any infalling 
subhalo remains intact, and we assume that this 
subhalo harbors a galaxy that will survive forever. 
This is tantamount to assuming that galaxies form 
in all sufficiently-large peaks in the primordial den- 
sity field and survive until today. 

Fric. Only - a model that only considers the effects 
of halo merging and dynamical friction. Subha- 
los never lose mass and can only be destroyed by 
sinking to the very centers of their hosts. 

Strip. Only - a model that only considers halo merging 
and mass loss and assumes no dynamical friction or 
central merging. Subhalos can lose mass and drop 
out of a mass threshold sample, but they cannot 
lose orbital energy and sink to the center of the 
host potential. 

Full - our full model treating halo merging, dynamical 
friction, and mass loss. This is the model that was 
developed in Z05 and validated against TV-body 
simulations. 

We run our models for host masses 3 in the range 
from log(Mhost//i _1 -M©) = 11.0 to 15.0 in steps of 
A(logMbost) = 0-1- F° r each of these masses, we run 
1000 statistical model realizations representing different 
realizations of the local density field and different halo 
merger histories. In this way, we sample the statistical 
properties of subhalo populations over the entire range of 
host halo masses relevant to galaxy-galaxy correlations. 
We repeat this process for host masses at z — as well 
as two past redshifts, z — 3 and z = 1, and two future 
redshifts, z = —0.6, and z = —0.9. 

4. EFFECTS OF SUBHALO DYNAMICS ON THE GALAXY 
CORRELATION FUNCTION 

4.1. Halo Occupation Distribution Statistics 

3 We note that wc use the "virial" definition of a halo in which a 
halo is defined as a spherical region of mean density equal to A v j r 
times the mean background density. For our cosmological model, 
A v j r = 337 at z = and approaches 178 at high z. 
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The galaxy correlation function may be consid- 
ered primarily a functio n of the galaxy HOD (e.g., 
iBerlind fc Wein berg 2002|). The prevailing cosmological 
model is now stringently constrained and may be con- 
sidered fixed for our purposes. Moreover, theoretical 
predictions of the abundances, clustering, and structures 
of host dark matter halos in the concordance cosmology 
are now well established. Consequently, we focus on the 
properties of the HOD and the manner in which the HOD 
determines galaxy clustering. 

We expect that each host halo of sufficient size con- 
tains one dominant, central galaxy associated with the 
host itself, as well as additional satellite galaxies that are 
associated with relatively large subhalos. Thus, the HOD 
of galaxies should resemble the HOD of all halos (hosts 
plus their subhalos), and suc h a model is bolstered by 
significant empirical support (iKravtsoy fc Klypinl IT999I; 
Colin et al.lll999tlKravtsov et alJl2004aHTasitsiomi et al.l 
2004HConrov et alJl200fih 7 As a result, we concentrate on 
the insight that can be gleaned about the development 
of the HOD of all halos, paying particular attention to 
the separate effects of halo mergers, dynamical friction, 
and mass loss. 

The left column of Figure [5] shows the mean occupa- 
tion number of host halos and subhalos as a function of 
host halo mass (N)m> &t z = 0. The three panels give 
results for halo samples defined by different mass thresh- 
olds. In the interest of simplicity, we assume that all host 
halos and surviving subhalos with masses M > M m j n 
harbor an observable galaxy. This assignment is simpler 
than those supported by detailed comparisons to data, 
which typically assume that all host halos and surviv- 
ing subhalos with masses M > M m i n (or some maxi- 
mum circular velocity) at the epoch of accretion harbor 
an observable galaxy. We proceed in this manner be- 
cause the subtleties discussed in the aforementioned liter- 
ature do not influence our primary points and may serve 
to obscure them. This is primarily because any mass 
threshold chosen at the epoch of accretion will have a 
second "destruction" threshold due to the finite resolu- 
tion of a given N-body simulation. This can alter clus- 
tering measurements, and since the aim of this paper 
is to present the qualitative trends responsible for the 
low-redshift correlation function, we use the simpler fi- 
nal mass approximation. We have confirmed that us- 
ing mass at accretion with our model reproduces the 
same general results. In a forthcoming paper we con- 
sider more sophisticated models to compare with data. 
In the top, middle, and bottom panels we show samples 
with log(M Bl i n //i- 1 M ) = 11.4, 11.7, and 12.3, respec- 
tively. These particular mass thresholds result in average 
galaxy number densities (see Eq. [3] ) equal to those in ob- 
served SDSS samples with r-band l uminosity threshold s 
of M r < -18.5, -19.5, and -20.5 (|Zehavi et al]l2005h . 
The four curves in each panel represent the four model 
modes described in § [3l and we calculate each curve from 
the mean of the 1000 model realizations. 

First, the black dot-dashed curves represent the No 
Effects model. As explained in § this model assumes 
that any halo that merges into a larger host system (and 
becomes a subhalo of that system) is thereafter unal- 
tered by dynamical effects in the host halo environment. 
Physically, this corresponds to the simple (and observa- 
tionally untenable) assumption that each subhalo above 



-^min brings with it an observable galaxy upon merging 
into the host and that this galaxy is not destroyed or 
dimmed by dynamical evolution within the host halo. In 
effect, each local peak of sufficient mass in the smoothed 
density distribution forms a galaxy and the galaxy can 
not be destroyed. Of course, we expect the mean halo 
occupation in this model to be high for all host masses 
as compared to the other models. 

Next, we turn to the curves depicting the individual 
effects of dynamical friction (Fric. Only, green dashed 
curves) and mass loss (Strip. Only, blue dotted curves). 
These dynamical mechanisms can destroy subhalos, but 
they cannot affect the host halo or central galaxy. This 
is why all the curves converge to the value (N)m = 1 
at low host masses. As a convenient shorthand, we refer 
to any subhalo that fell into its host system with a mass 
M su i> > M m ; n , but then merged with the central host 
galaxy or lost sufficient mass to fall below this threshold, 
as destroyed. This does not mean that the subhalo has 
become unbound, but merely that it has either merged or 
no longer has a bound mass above some minimum mass 
threshold. 

Dynamical friction acting alone destroys subhalos by 
causing them to sink to the centers of their hosts and 
"merge" with it. This mechanism alone causes a 20- 
35% decrease in the mean number of surviving satel- 
lites for all host masses as compared to the No Effects 
model. The fractional decrease in subhalos depends only 
weakly on host mass, but a comparison of the differ- 
ent mass threshold panels shows a modest dependence 
on subhalo mass, with smaller mass subhalos being de- 
pleted more. These trends are counter-intuitive because 
the dynamical friction force is an increasing function of 
-Msub /Mhosti the mass ratio bet ween the subhalo and its 
host (|Binnev fc Tre maiiic 2008). One might expect the 
depletion of subhalos to be larger for smaller host masses 
at fixed subhalo mass or for larger subhalo masses at fixed 
host mass. Our results differ from this expected behav- 
ior for two reasons. First and foremost, low-mass-ratio 
mergers tend to occur at higher redshifts than high-mass- 
ratio mergers. At higher redshifts, host halos are signifi- 
cantly smaller than they are at present, so high-redshift 
mergers probe only the dense interiors of contemporary 
host halos and evolve approximately according to the 
subhalo-host mass ratio at the redshift of the merger. 
These early-merging subhalos also have a longer period 
of time during which to evolve. Second, our models in- 
clude subhalos-of-subhalos. As we move to larger host 
masses at fixed subhalo mass or smaller subhalo masses 
at fixed host mass, more subhalos are subs-of-subs that 
have much higher mass ratios with their immediate hosts. 
These effects result in the trends we see in Figure [5] 

Mass loss is significantly more effective at "erasing" 
subhalos than dynamical friction. Mass loss processes 
can effectively "destroy" subhalos because many lose suf- 
ficient mass to fall below the threshold of a sample. This 
mechanism typically drives an 80-85% decrease in the 
number of objects above a given mass threshold com- 
pared to the No Effects model. Again, the fractional 
decrease in subhalos is nearly independent of host mass, 
but it shows a slight dependence on subhalo mass, with 
smaller mass subhalos being destroyed more efficiently. 

Finally, the Full model (red, solid curves) includes the 
effects of both subhalo mass loss and orbital decay by 
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Fig. 2. — Left Panels: Mean number of all halos (hosts plus subhalos) predicted by our subhalo model as a function of host halo mass, at 
redshift 2 = 0. The three panels show results for three mass threshold values: log(M m i n /fc _1 M ) = 11.4, 11.7, and 12.3. The four curves 
in each panel correspond to the four models described in § [3] No Effects considers no gravitational effects on subhalos as they orbit inside 
their host halos (black dot-dashed curve); Fric. Only considers only the effects of dynamical friction (green dashed curve); Strip. Only 
considers only the effects of mass loss (blue dotted curve); Full considers both dynamical friction and mass loss (solid red curve). Right 
Panels: The correlation function of all halos predicted by our subhalo model. f(r) is computed from the halo model using the occupation 
statistics shown in the left panels. The figure shows that dynamical effects (especially mass loss) are needed in order to reduce the number 
of subhalos sufficiently and produce a power-law correlation function. 



dynamical friction. These processes do not simply sum 
together. As a subhalo sinks deeper into its host po- 
tential well due to dynamical friction, it experiences a 
stronger tidal field and is thus more efficiently stripped 
of its mass. Conversely, less massive subhalos are less 
susceptible to orbital decay via dynamical friction. A 
comparison of the Full model to the Strip. Only model 
shows that including dynamical friction causes an addi- 
tional ~ 15% depletion of substructure. Mass loss is by 
far the dominant cause of subhalo destruction. Overall, 
Figure [2] shows that dynamical effects reduce the number 
of subhalos by ~ 90% compared to the number of dis- 
tinct mergers that occur during the formation of a host 
halo. 

4.2. Constructing the Correlation Function 

We use the halo model outlined in § 12.11 to com- 
pute the correlation function pre dicted by our subhalo 
model. Specifically, we u se the Uenkins et al.l (|2001l ) 
mass function and we follow iTinker et all (|2005l) m using 
the ISmith et al.l (|2003t ) form ula for the non-linear mat- 
ter power spectrum and the ITinker et al.l (|2005h scale- 
dependent halo bias relative to the non-linear power spec- 
trum. We derive HOD statistics from our subhalo mod- 
els as exemplified by the previous section, and we com- 
pute the pair separation distributions by assuming that 
the radial distributions of satellites follows an NFW pro- 
file for simplicity. In actuality, the subhalo distributions 
in both our models and TV-body simulations are slightly 
shallower than NFW (see Z05 for model and simulation 
results). We adopt the NFW profile for analytical con- 
venience as deviations from NFW are small and only 
influence correlation functions notably on scales signifi- 
cantly smaller than r ~ 100/i -1 kpc (e.g., Z05; also see 
IWatson et al.| [2010 for a demonstration of this point re- 
garding satellite galaxies). 

The right column of Figure [5] shows the host+subhalo 
correlation functions computed in this manner from the 
HODs predicted by our subhalo model. In the No Effects 
case, where no subhalos are destroyed, £(r) is very dif- 
ferent from a power law, having a one-halo term that is 
too large relative to the two-halo term so that a distinct 
feature is present in £(r) on scales r ~ 2 Mpc. In fact, 
comparing this to Figure [IJ we see that it is very similar 
to the dark matter correlation function. This is perhaps 
not surprising because subhalos in this model behave as 
massive test particles that cannot be altered. As dynam- 
ical effects are included and substructure is consequently 
depleted, £(r) drops on all scales. Recall that only subha- 
los (hosting satellite galaxies) can be destroyed and the 
number of host galaxies remains unaltered. The fraction 
of all objects that are satellites therefore decreases. As 
we discussed in § 12. H for a fixed population of central 
galaxies, the one-halo term drops in approximate pro- 
portion to the number of satellite galaxies. So as the 
number of satellites declines, so does the number of pairs 
within halos relative to the total number of pairs and the 
one-halo term declines. 

Large-scale clustering is less sensitive to changes in the 
satellite galaxy population. The two-halo term drops 
because subhalos tend to populate more massive hosts 
(as in the left column of Figure [2]), so the average host 
halo mass of a sample decreases as subhalos are depleted. 
The large-scale clustering strength of halos increases with 



halo mass, so this depletion results in weaker large-scale 
clustering. The variability of the two-halo term is rela- 
tively mild because the halo bias is not a ra pidly- varyin; 



pidiy-varying 
(jTinker et all 



funct ion of halo mass near M 
l2005h . 

With enough depletion of substructure, the one- and 
two-halo terms align and result in a nearly power-law 
shape. This is exactly what happens in Figure [2j In our 
Full subhalo model the correlation function is roughly a 
power law. To obtain a nearly power-law galaxy corre- 
lation function, it is necessary that a majority of early 
galaxies and proto-galaxies that merge to form a massive 
system at low redshift be destroyed through either cen- 
tral mergers or mass loss. Our results suggest that mass 
loss is mainly responsible for this depletion, while dynam- 
ical friction and central galaxy mergers play a compara- 
bly small, supporting role. Incidentally, this picture im- 
plies that infalling satellite galaxies lose significant stel- 
lar mass so that they provide an important source of the 
diffuse intracluster light observed in galaxy groups and 
clusters and this picture i s consistent with observations 
(|Purcell et al.l l20Q7l I2Q08D . Comparing our correlation 
function results for the three different mass thresholds, 
we note that £(r) is closer to a power law for lower-mass 
samples. We revisit this point in the next section. 

We now return to the mean occupation statistics shown 
in the left panels of Figure [2j The so-called "plateau" re- 
gion of the HOD is the flat region at (N)m = 1, where 
host halos are more massive than M m ; n , but not yet mas- 
sive enough to host subhalos above our mass threshold 
4 . As substructure is depleted, the prominence of this 
plateau increases. The "length" of the plateau in the 
HOD can be expressed as the ratio between the mass of 
a halo that hosts a single satellite on average, Mi, to the 
minimum mass req uired to host a central galaxy, M m { n . 
IZehavi et al.l ([2005) fit an HOD model to the measured 
correlation function of SDSS galaxies and found that a 
ratio Mi/M m j n ~ 23 is consistent with clustering data, 
nearly independent of galaxy luminosity. In other words, 
a consistent picture is one in which the entire HOD shifts 
to higher masses in a self-similar manner, with Mi/M m j n 
fixed, in order to accommodate hi g her-lum inosity sam- 
ples. Remarkably, iKravtsov et al.1 (j2004aT l studied this 
for subhalos in a high-resolution iV-body simulation and 
found that Mi/M m j n ~ 2 0, reg ardless of M m i n as well. 
Meanwhile, ITinker et all ([2005) fit a slightly more com- 
plex HOD model to the SDSS data and found that 
Mi /Mmin ^25 for galaxy samples with luminosities less 
than L*, but decreases to Mi/M m j n < 5 to accommodate 
the highest-luminosity samples (absolute r -band magni- 
tudes M r < —21). The new analysis bv IZehavi et al.1 
(|2010f ) also finds this trend with M 1 /M min ~ 17 for 
M r > —20.5 and much lower values for higher luminosity 
galaxies. For the purpose of comparison, the Full sub- 
halo model shown in Figure [2] predicts Mi/M m i n ~ 40 
for the low-mass samples of log(M m i n //i -1 M©) = 11.4 
and log(M min /^- 1 M Q ) = 11.7, and Mi/M min ~ 30 for 
the higher-mass sample of log(M min / 'h~ x M©) = 12.3. 

These results suggest that getting the length of the 

4 Roughly speaking, the most massive subhalo within any host 
is a few percent of the mass of the host halo (e.g., Z05). This is 
the case with, for exam ple, the Large Magellanic Cloud within the 
halo of the Milky Way l IBusha et aUl20iafl . 
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HOD plateau right may be a key ingredient needed to 
establish a power-law correlation function and this has 
been part of the interpretation in the literature. The 
importance of Mi/M m - m stems from the fact that most 
one-halo pairs reside in halos with average satellite num- 
bers (N s )m J$ 1, so modeling the HOD at relatively low 
satell ite occupation numbers is critical (see lConrov et all 
|2006() . We investigate this further in § [6] 

5. MASS AND REDSHIFT DEPENDENCE OF THE 
CORRELATION FUNCTION 

5.1. Dependence on Mass 

While dynamical processes act in a manner to deplete 
substructure and push £(r) toward a power law at all 
mass thresholds, it is evident that deviations from a 
power law are stronger with increasing host mass. Fig- 
ure [3] shows the correlation functions predicted by our 
Full subhalo model for four different mass thresholds, 
ranging from log(Mmin/ h~ x Mo) = 13.5, correspond- 
ing to bright galaxies such as Luminous Red Galaxies 
(LRGs), down to log(Af min //i _1 M Q ) = 10.5, correspond- 
ing to dwarf galaxies. While the correlation function 
of the "dwarf" sample is a near power law, that of the 
"LRG" sample exhibits strong departures from power- 
law behavi or. This trend has been detected with SDSS 
galaxies bv lZehavi et al.l (|2005h who found evidence that 
a power-law model provides a better fit to low-luminosity 
galaxies than high-luminosity galaxies. Halo and subhalo 
clustering exhibits the same trend. More massive halos 
contain slightly more of their bound masses in substruc- 
ture relative to less massive halos, but this is a compa- 
rably small effect (Z05) and drives only ~ 30% of the 
mass-dependence of the one-halo term in Figure [3] At 
fixed redshift, the departure from a power-law at high 
mass (high luminosity) is caused by the relative rareness 
of high- mass host halos (see § l2.1l for interpretive discus- 
sion) . 

5.2. Dependence on Redshift 

Substructure abundances vary with time. Infall of new 
subhalos acts as a "source" of halo substructure. The 
rate of mergers of halos into larger systems is a function 
of redshift that typically peaks at redshifts z ~ 1 — 3 
in th e halo mass ra nge of interest and declines thereafter 
(Z05: lZentner1l2007f) . Once a subhalo merges into a larger 
host halo, dynamical friction shrinks its orbit and the 
subhalo loses mass. Given enough time, the subhalo will 
eventually lose enough mass to fall below M min or merge 
with the central galaxy and lose its identity. The balance 
between the halo merger rate and the rates of destruc- 
tive processes (which occur on a halo dynamical time) 
determine the redshift dependence of halo substructure. 

Figure [4] shows the redshift evolution of the mean halo 
occupation number and resulting correlation functions. 
The layout of Figure@]is similar to that of Figure^ with 
(N)m shown in the left panels, £(r) shown in the right 
panels. However, in Figure 0] all results are for the Full 
subhalo model, and the various lines denote quantities 
evaluated at different redshifts, z = 3, 1, 0, —0.6, —0.9 
(where negative redshifts correspond to future epochs). 
Moreover, in each panel the correlation functions are 
scaled by a power law to better highlight departures from 
a power-law shape. 
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Fig. 3. — Correlation function of all halos (hosts plus subhalos) 
predicted by our subhalo model at redshift z = 0. The four curves 
show £(r) for four mass threshold samples and the threshold values 
-Mmin (in units of h~ Mq) are listed in the panel. The figure 
shows that f(r) breaks more and more from a power law for higher 
mass halo samples, which correspond to higher luminosity galaxy 
samples. 



The left-hand panels of Figure [4] show that the average 
number of subhalos within hosts of a given mass starts 
out high at early times and begins to decrease after z = 3, 
as merger rates decline. By the present epoch (z = 0), 
the number of subhalos has dropped by ~ 25 — 30% rel- 
ative to what it was at z = 3. One Hubble time into the 
future [z = —0.6), the abundance of substructure has 
dropped by ~ 60%. This is because the rate of merging 
as a source for new subhalos declines rapidly. This de- 
crease in the merger rate is dictated in large part by the 
quenc hing of structure g rowth by the cosmological con- 
stant (|Carroil et al.|[l992 ). but also because most halos of 
interest are below the typical collapsing mass, which ap - 
proaches ~ 10 14 h^M® in the future (|Zentnerll2007l) . 
Meanwhile, destructive processes continue to operate on 
orbiting halo substructure for several additional dynami- 
cal times. Three Hubble times into the future (z = —0.9) 
the average halo occupation has dropped by ~ 90%. As 
with our previous results, the fractional decrease in sub- 
halo abundance appears to be roughly independent of 
host mass, meaning that the slope of the HOD in the 
high-(iV)jvf limit is not significantly altered by evolution. 
The amplitude of (N)m declines considerably, resulting 
in increasing Mi/M m i n , or a "lengthening" of the HOD 
plateau with time. This behavior is strikingly similar to 
that seen in Figure[2]in the sense that turning on dynam- 
ical effects at a fixed redshift has a qualitatively similar 
impact as evolving forward in time, and the effects on 
the correlation function are similar. 
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Fig. 4. — Left panels: Mean number of all halos (hosts plus subhalos) predicted by our Full subhalo model as a function of host halo 
mass, at five different redshifts. The three panels show results for three mass threshold values: log(M min //i _1 M ) = 11.4, 11.7, and 12.3. 
The five curves in each panel correspond to the redshifts z = 3,1,0,-0.6,-0.9 (negative redshifts correspond to future epochs). Right 
panels: Correlation functions corresponding to the halo samples shown in the left panels. In each case, £(r) has been scaled by a power law 
in order to clearly show departures from a power-law shape. The figure shows that the number of subhalos steadily decreases from high 
to low redshift, causing the correlation function to evolve from not being a power law at high redshift, towards having a nearly power-law 
shape at the present epoch, and once again deviating from a power law at future epochs. 
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Fig. 5. — The correlation function of all halos (hosts plus subha- 
los) predicted by our Full subhalo model as a function of redshift, 
for a single mass threshold sample log(M m ; n //i — 1 Mq) = 12.3. 
Each panel shows § (r ) for a different redshift (solid curve) , as well 
as the one-halo (dashed curve) and two-halo (dotted curve) terms. 
The figure shows that the one-halo term evolves strongly with red- 
shift and only at z = strikes the right balance with the two-halo 
term to result in a power law. 



Turning to the right panels, £(r) is shown at each red- 
shift scaled by an r -17 power law in order to emphasize 
features in the correlation function. Starting at z = 3 
(long dashed curves), £(r) is very far from a power-law, 
with a slope that is much steeper on small scales. At 
z = 1 (short dashed curves) the break from a power 
law is less pronounced, but it is still significant. These 
results are qualitatively co nsistent with clustering mea- 
sure ments at high red shifts ()Coil et al.ll2006t lOuchi et all 
120051 iLee et al.ll2006h . At z = (solid curves), the cor- 
relation function is approximately a power law, though 
there is still a mild, discernible feature at the transition 
scale between the one- and two-halo terms. In the future, 
£(r) once again breaks from a power law. At z = —0.6 
(dot-dashed curves), departures from a power-law shape 
are about as strong as they were at z = 1 . Three Hubble 
times into the future, at z — —0.9 (dotted curves), the 
departures from a power law are significant and repre- 
sent a dramatic reduction in the relative contribution of 
the one-halo term. 

Figure [5] focuses on the log(M min //i~ 1 Af ) = 12.3 
threshold sample and shows the correlation function at 
four different redshifts, while also showing the one-halo 
and two-halo terms explicitly. Figure [5] clearly demon- 
strates how a delicate balance is needed between the two 
terms in order for £(r) to achieve a power-law shape. The 
two-halo term exhibits modest variations from panel to 
panel, with a range of about a factor of ~ 3. The de- 
creased large-scale clustering at z > is due to the linear 
growth of perturbations with time, but this is always kept 



modest because t he increasing b ias of halos of fixed mass 
with redshift (see IZentnerl 12007) compensates for large- 
scale structure growth. At z < 0, the slight decrease in 
two-halo clustering is due to the decay of halo bias once 
halo growth slows (|Frylll996T ). 

The variation in the one-halo term is significantly 
larger, as our earlier discussions suggest, and changes by 
a factor of ~ 45 — 150 (depending on scale), equivalent 
to ~ 15 — 50 times the variation in the two-halo term. 
At high redshift, the relative rareness of host halos and 
the large amount of substructure cause £ (r) to be boosted 
significantly in the one-halo regime as shown in the z = 3 
panel of Figure El At z = 0, just the right amount of 
substructure has been depleted to strike a near balance 
between the one-halo and two-halo contributions. In the 
future, the continual destruction of subhalos suppresses 
the one-halo term, driving £(r) away from a power law 
again. By z = —0.9, the depression in small-scale clus- 
tering is striking. 

Some of the evolution of £(r) on small-scales comes 
from the fact that halos large enough to host luminous 
galaxies become increasingly rare as redshift increases. 
The characteristic collapsing mass is a rapidly decreas- 
ing function of redshift and is only M* « 10 9 /i _1 M 
at z = 3. In the relevant regime, the strength of the 
one-halo term grows in approximate proportion to the 
number of satellite galaxies and in inverse proportion to 
the number of host halos of appropriate size (see § 12. ip . 
so the relative paucity of host halos at high redshift also 
drives strong one-halo clustering because Fig.[5]describes 
samples of fixed absolute mass threshold. However, it is 
subhalo abundance that has the larger influence on the 
redshift dependence of clustering. We have computed 
the correlations of Figure [5] using samples in which M m j n 
varies with redshift so as to maintain a constant number 
density of halos. These samples are less subject to the 
gross evolution of the halo mass function. We find all of 
the same qualitative results for this case, though the two- 
halo term varies by a factor of ~ 4, while the variation 
in the one-halo term is limited to a factor of ~ 12 — 80 
(again, depending on scale), resulting in a variation in the 
one-halo term that is ~ 3 — 20 times larger than that of 
the two-halo term. Moreover, we have re-computed cor- 
relation functions using a combination of the predicted 
low-redshift HODs alongside the high-redshift mass func- 
tions in order to isolate the contribution due to the mass 
function and HOD evolution. The majority the redshift 
dependence of £(r) on small scales is due to the evo- 
lution of subhalo abundance. To maintain a power-law 
correlation function at high-redshift would require fewer 
subhalos per host than at z = in order to compensate 
for the relative rareness of host halos at high-redshift. 
In fact, hosts at high redshift have a larger number of 
subhalos of any given mass so these effects reinforce one 
another, leading to a strong deviation from a power-law 
£(r) at high redshift. 

We have already described the reasons that the one- 
halo and two-halo terms behave so differently under 
changes in the HOD. To reiterate, on large scales, £(r) 
is essentially a weighted average of the clustering of host 
halos, where {N)m provides the weighting (see the in- 
tegral in Eqs. [6] and [7], note that X(k,M) « 1 for 
k < l/i? V j r ). The possible variability in £(r) on large 
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scales is limited because it is always bound by the lim- 
ited variation in the clustering of host halos. As we dis- 
cussed in § 12.11 the difference in the large-scale bias of the 
largest relat ive to the smallest halos is at most a factor 
of ~ 3 (e.g.. lTinker et aDl2005f ). Significant variations in 
large-scale clustering require dramatic variations in the 
HOD at high mass, which are not expected on theoret- 
ical grounds and are not mandated by data. However, 
on scales smaller than the size of individual host halos, 
£(r) can vary dramatically depending on the HOD. For 
example, in the extreme case of only one object per host 
halo, there will be zero pairs within halos and the one- 
halo term will vanish. For large numbers of satellites, the 
one-halo term will be significantly larger than a power- 
law extrapolation of the two-halo term to small scales. 

The sensitivity of the one- halo term to the HOD, cou- 
pled with the relative insensitivity of the two-halo term, 
means that achieving a power-law correlation function re- 
quires fine-tuning in the number of satellite galaxies per 
halo. The satellite galaxy abundance naturally evolves 
with redshift, so £(r) can only be a power law during 
those epochs when substructure has evolved to join the 
one-halo term to the two-halo term. Of course, it may 
be possible for features in the host halo mass function or 
bias relations to conspire to compensate for substructure 
evolution, but such features would somehow need to be 
coordinated with low-redshift structure growth. A dif- 
ferent way to state this is that the halo mass function 
and halo bias depend on the statistics of the linear den- 
sity field, and do not "know" about the non-linear galaxy 
formation and gravitational processes that occur within 
halos. It would be quite strange if their evolution were 
somehow connected with the evolution of satellite galax- 
ies in virialized hosts. It seems to be a coincidence that 
the epoch of near power-law clustering of typical galaxies 
lies near z = 0. 

5.3. The Balance Between Accretion and Destruction 

We have just seen how the depletion of substructure 
over time leads to evolution in the correlation function 
such that it becomes a power law at the present epoch. 
However, what drives substructure depletion? We expect 
that most subhalos will lose significant amounts of mass 
or merge with the central galaxy given sufficient time. 
However, this will be compensated to some degree by 
the infall of new subhalos. If the rate at which satellites 
are accreted is greater than the rate at which they are 
destroyed, then the net amount of substructure will grow 
with time. The evolution in the number of subhalos (and 
hence the correlation function) depends on the balance 
between accretion and destruction. Z05 give a related 
discussion of accretion and destruction in their § 4.4 and 
the perspective we adopt here complements Z05. 

In Figure [6] we illustrate the competition between ac- 
cretion and destruction in host halos of mass M = 
10 134 /i _1 M . To measure the accretion rate (dashed 
curve), we count all subhalos with masses greater than 
lO 1:L ft. _1 M0 that accrete onto these hosts in finite time 
intervals. For the destruction rate (solid curve), we 
count the number of these same subhalos that drop be- 
low 10 11 h~ x Mq during the time intervals. The accretion 
rate minus the destruction rate will then give us the net 
rate of change in the number of subhalos per unit time. 

Figure[6]shows that the accretion rate quickly grew and 
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Fig. 6. — The accretion versus destruction rate of subhalos over 
cosmic time, as predicted by our full subhalo model. The accretion 
rate shown is the number of subhalos per Gyr that merge into a 
host halo of mass log(M/h~ 1 M < 7 ) ) = 13.4. The destruction rate 
is the number of these same subhalos per Gyr that are destroyed 
(i.e., their mass drops below some threshold value). The two rates 
equalized when the Universe was ~ 6Gyr old (at z ~ 1). Before 
z = 1, the net number of subhalos increased with time, whereas at 
later times the net number decreased with time. The figure shows 
how the balance between accretion and destruction changes with 
redshift, which explains why the correlation function can only be 
a power law at a single epoch. 

reached a peak at z ~ 2— 3. Since this peak, the accretion 
rate has been steadily declining and is close to zero at the 
present epoch. The decline in merger rate s is partly due 
to the shape of the power spec t rum (see lLacev fc Cold 
Il993t ISomerville fc KoIattlll999t IZentnerll2007ft . but the 
driving force for the recent fast decline in the merger rate 
of halos is the reduced rate of structure growth caused by 
accelerated cosmic expansion. The destruction rate also 
peaked at z ~ 2 — 3 and lags the accretion rate because 
most destruction happens over a period of several dy- 
namical times. Figure |5] clearly shows that the accretion 
rate has been dropping faster than the destruction rate 
since their peaks, with accretion an d destruction roughl y 
balancing just below z ~ 1 (see also lStewart et aLi r2009'). 
This means that the number of subhalos in hosts that 
grow to a mass of M = 10 13 ' 4 ft, M Q by z = increased 
until z ss 1 and has been declining ever since, despite the 
fact that the virial masses of these halos have been grow- 
ing. The general trend toward reduced substructure at 
low redshift explains the behavior exhibited in Figure 2J 
The correlation function is close to a power law at the 
present epoch because the balance between accretion and 
destruction over time has led to the requisite abundance 
of substructure today. 

6. ACHIEVING A POWER-LAW CORRELATION 
FUNCTION 

We now step back from making predictions using our 
specific subhalo model and undertake a general explo- 
ration of the properties of the HOD that yield nearly 
power-law correlation functions at different masses and 
rcdshifts. The HOD characterizes the number and spa- 
tial distribution of galaxies within dark matter halos. It 
is typically specified with a handful of parameters that 
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are c onstrained using galaxy clustering measurements 
(e.g.. iMagliocchetti fc Porcianiir2 003: Zc havi et al.1 120051 : 
iTinker et al.1 120051: IZheng et al.1 l2007t) . We choose an 
HOD model that is motivated by theoretical predictions 
from hydrodynamic simulations, semi-an alytic models, 
and high-resolution N-body simulations (IBerlind et al.1 
120031: IKravtsov et alj|2004at IZheng et al.ll2005D . Accord- 
ing to this model, halos above some threshold mass con- 
tain a single "central" galaxy plus a number of "satel- 
lite" galaxies. The number of satellites in any given halo 
is drawn from a Poisson distribution whose mean is a 
power- law function of host halo mass. The central galaxy 
is placed at the center of the host halo, while the satellites 
are spatially distributed according to an NFW density 
profile. Specifically, we adopt an H OD paramete r izatio n 
that is similar to the one used bv ITinker et~aH ([2005j) . 
This is a simple, yet powerful model in which the num- 
ber of central galaxies is modeled as a step function, 



1 if M > M min 
if M < M min 



(8) 



while the mean number of satellites follows a power-law 
with an exponential cutoff at low mass, 



{N sat )M = TT exp 




-M, 



M 



(9) 



The parameters in the model are as follows. 



1. M m i n is the minimum host halo mass to contain a 
central galaxy. 

2. Mq is the host halo mass below which satellite 
galaxies are exponentially supressed. 

3. Mi is the host halo mass to contain, on average, 
one satellite galaxy. 

4. a is the index of the power-law relation between the 
mean number of satellite galaxies and halo mass. 

Previous studies have shown that the power-law 
index a ps 1 for subhalos and simulated galaxies 
(IKravtsov et al.ll2004at IZheng et alJl2005l IZentner et al.1 
I2005D. as well as observed galaxies dimmer th an L* 
(jZehavi et aI1l2005D . leading ITinker et al.1 (|2005l ) to set 
a = 1 throughout their analysis. However, we allow 
a to vary because the correlation function is sensitive 
to it and, although it may be near unity when model- 
ing observed data, it may need to deviate from unity to 
yield a power-law correlation function at high redshifts. 
On the other hand, £(r) is not sensitive to Mq, conse- 
quent ly we fix its value by adopting the IConrov et al.1 
(|200ll M - M x relation, 

log(Mo//i -1 M Q ) = 0.76 log(Mi/h~ 1 M (3 ) + 2.3. (10) 

The result is an HOD model with only three free param- 
eters: M min , Mi, and a. 

The two-halo term of £(r) depends on the mean oc- 
cupation (N)m — (N ccn + A^t )m> which is equal to 
1 + (A 7 sa t)M f or M > M min . The one-halo term also 
requires the second moment of the occupation distri- 
bution (A sat (A r sat — 1))j\/i so characterizing the mean 
occupation is not sufficient. We assume that the num- 
ber of satellites follows a Poisson distribution, for which 



(N sat (N sat - 1))m = (N sa t} 2 M . Our Full model devi- 
ates mildly from a pure Poisson d istribution (see Fig. 7 
of Z0 5, and recent simulations of iBovlan-Kolchin et al.1 
l201Ct that find similar deviations from a Poisson distri- 
bution), but the effect of this deviation on £(r) is minor 
(Fig. 16 of Z05). We also note that there are any num- 
ber of possible parametrizations for {N sat )M to choose 
from besides the one adopted here. We have found that 
mildly different parametrizations that exhibit the same 
basic features and are c onsistent with conte mporary data 
(e.g., the one used bv IZehavi et al.l 120051 ) yield similar 
conclusions. 

We consider the HOD parameter space that yields a 
power-law correlation function for three galaxy samples 
of fixed number density n g , at three different redshifts 
z = 0, 1, 3. Fixing number density is a way to compare 
similar samples at different redshifts because the high- 
redshift sample is more likely to represent the progeni- 
tors of the low-redshift sample than it would in the case 
of mass threshold samples. We choose number densities 
equal to n g = 0.02, 0.01, and 0.003 /i 3 Mpc" 3 , which cor- 
respond to three z ~ 0, volume-limited, r-band threshold 
samples in the SP SS: M r < -18.5, -19.5, and -20.5 
(jZehavi et al.ll2005| ). 

For a given number density and redshift combination 
(e.g., n g = 0.02 h 3 Mpc~ 3 at z = 1), we create a 50 x 50 
grid of Mj — a parameter combinations. For each pair 
of Mi and a on this grid, we use Equation ([3]) to find 
the value of M m - ln that is needed to enforce the desired 
number density. In this manner, the 2,500 HOD mod- 
els on the grid represent galaxy samples with the same 
number densities, but different HODs. We then com- 
pute the first and second moments of the mean galaxy 
occupation using equations <j8j) and © , and use the halo 
model described in § 2.1 to construct We assign 10% 
errors on all scales to £(r), as such errors are roughly 
consistent with jackknife re-sampling errors in current 
clustering measurements (|Zehavi et al.ll2005h . and we fit 
a power-law function to all 2,500 correlation functions. 
We perform our fits using a Markov Chain Monte Carlo 
(MCMC) analysis in which we vary the slope and corre- 
lation length of the fitted power-law. We then find the 
minimum \ 2 value for a power-law fit to £(r) for any 
given Mi — a combination. This allows us to approxi- 
mate the HOD parameter space in which £(r) is consis- 
tent with a power law at a level similar to contemporary 
observations. For two free parameters, the 68.3% (ler), 
95% (2a), and 99.6% (3a) likelihood regions correspond 
to values of reduced x 2 < 1.15, 1.61, and 2.06, respec- 
tively. 

Figure [7] shows the contours generated from the afore- 
mentioned procedure. Each row in the figure represents 
a different n g value and each column corresponds to a 
different redshift. The "satellite fraction" (/ sa t, the frac- 
tion of all galaxies that are satellites, see Eq. [5]) is 
relevant to the shape of the galaxy two-point correla- 
tion function. Therefore, over-plotted in each panel are 
curves of constant / sat (the labeled, solid, black curves). 
To compare these results with measurements from ob- 
served galaxies, in each panel we also show best-fit Mi 
and a values from published halo model fits to mea- 
surements of £(r) using galaxy samples with the same 
number densities at the same redshifts. Squares and 
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Fig. 7. — Exploration of the HOD parameter space that yields a power-law £(r), as a function of redshift and sample number density. 
Each column of panels shows results for a different redshift [z = 0, 1,3). Each row of panels shows results for a different sample number 
density (n g = 0.02, 0.01, 0.003/i 3 Mpc~ 3 ). We adopt the four-parameter HOD model shown in equations [8] [9] and 1101 Each panel shows 
the parameter space probed by a, the slope of the mean occupation number of satellites, and Mi , the halo mass that contains on average 
one satellite galaxy. For each pair of a and Mi values, we find the value of M m i n that yields the desired galaxy number density. We then 
use the halo model to compute §(r) for that set of HOD parameters. We do this on a 50 X 50 grid of a — Mi parameter combinations. 
We fit each correlation function to a power law, and the shaded contours represent the 68.3%, 95% and 99.6% power-law likelihood (green, 
blue, red contours). Also shown are contours of constant satellite fraction (solid black curves). The red X in each panel shows the HOD 
parameters predict ed by our Full sub halo model. For comparison, we also show results from HOD modeling of real ga laxy samples from 
the SDSS at z = l lZheng et alJ2007l T05 - magenta boxes and grey asterisks); and DEEP2 at z = 1 (Zheng ct al. 2007 - purple triangles). 
Finally, we show the simulation results of IConrov et al. (2006) that are designed to model SDSS, DEEP2, and Lyman-break galaxies at 
2 = 0, 1, and 3, respectively (cyan diamonds). 



triangles represent the best-fit parameter values from 
IZheng et all ([2001 who fit SDSS (z = 0) and DEEP2 
= 1) data (ZSDSS, ZDEEP2) and asterisks repre- 
sent the best-fit [Tinker et al.l (120051) values for SPSS data 
(T05). Diamonds represent the IConrov et al.l (|2006f ) val- 
ues for SDSS, DEEP2, and the z = 3 Subaru data of 
Lyman-break galaxies (C06). 

The best-fit parameter combinations should be re- 
garded as best-fit "regions" , because there are er- 
rors associated with the derived parameters (e.g., the 



IZheng et all (|2007f ) SDSS a and M x error s at each lu- 
minos ity are of order 10%). We note that I Tinker et al.1 
(2005) considered several possible values of a%, but we 
show their results for cr§ = 0.9 to be consistent with the 
cosmological model used in the other studies. Finally, in 
each panel we show the HOD parameters predicted by 
our Full subhalo model (marked by an "X" ) for samples 
with mass thresholds that yield the desired number den- 
sity. The Full model gives (-/V sat }M and we fit this with 
Equation ^ to obtain best-fit values of Mi and a. 
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Several interesting conclusions can be drawn from this 
figure. 

1. The region of HOD parameter space that yields 
a power-law £(r) drifts to lower values of both 
Mi and a with increasing number density. These 
trends increase the satellite fraction as number den- 
sity increases to compensate for the relative reduc- 
tion in the one-halo term compared to the two- 
halo term induced by moving to a lower-mass, more 
abundant halo sample. 

2. The values of a that result in the best power laws 
drift higher with increasing redshift in an effort to 
boost the two-halo term by placing galaxies in mas- 
sive, highly-biased halos. In general, it is difficult 
to arrange a power law at z > 3 for these three 
number densities. 

3. As might be expected from our previous discus- 
sions, there is a relatively narrow range of / sa t 
for the best-fit power-law space at each redshift. 
At z = 0, the space that is consistent with a 
power-law with 10% errors on the data lie near 
/sat ~ 0.1 — 0.15. The direct fits to observational 
data lie near / sat = 0.2 — 0.3. At z = 1, the 
power-law region is shifted to / sat ~ 0.05 — 0.1, 
while at z = 3 the power-law region is even lower, 
/sat ~ 0.01 — 0.02. Note that the power-law re- 
gions are not precisely aligned along constant- / sa t 
contours, particularly at low redshift and low num- 
ber density, indicating that other factors, such as 
host halo abundances and the physical sizes of host 
halos, contribute to the power-law nature of £(r). 
However, at high redshift and low number density, 
the power-law regions become more nearly co-linear 
with contours of constant / sa t over a range of a val- 
ues. 

4. The SDSS (z ~ 0) best-fit points lie near the power- 
law contours, but not within these likelihood re- 
gions. This is not surprising as the SDSS measure- 
ment is more precise over a wide range of scales 
than the ~ 10% errors we have assumed and the 
observed £(r) is now known to exhibit very small, 
but statistically-signific ant deviations from a power 
law ()Zehavi et al.ll2004D . 

5. As predicted from Figure El the fits to observa- 
tional data lie further from the power-law regions 
as we move to lower n g (higher luminosity) sam- 
ples. At fixed redshift, this is driven largely because 
the host halos of these galaxies become increasingly 
rare. However, it is worth noting that the growth 
of the one-halo term with increasing M m - ln is re- 
inforced by an increase in satellite abundance at 
fixed scaled mass M su b/Mh os t as Mh os t increases, 
accounting for ~ 30% of the rise. This increase 
satellite abundance with Mhost arises because more 
massive host halos assemble more recently, leaving 
less time for the evolution of substructure and less 
satellite destruction (Z05). The relative time avail- 
able for satellite evolution is an important part of 
determining the power-law nature of the correla- 
tion function. 



6. The fits to observational data lie near the power- 
law regions at z = 0, but grow more distinctly 
separated with increasing redshift. This evolution 
is driven by satellite fractions at high-z that are 
too large to be consistent with power-law cluster- 
ing. This supports our basic picture that satellite 
destruction over cosmic time is needed to achieve 
a power law £(r), and that the observed low- 
luminosity, low- redshift £(r) is a coincidence. 

7. The HOD values predicted by our Full subhalo 
model are similar to all of the observed data fits 
at all redshifts. This is a remarkable result con- 
sidering our model treats only subhalos and not 
galaxies explicitly. We explore more complicated 
associations of galaxies and subhalos in a follow-up 
study. 

8. Our subhalos, as well as all observational data, re- 
veal values of a ~ 1 for all reds hifts, in accord 
with p revious theoretical re sults (jKravtsov et al.1 
2004ai IZheng et all 120051 : IZentner et all 120051 : 
Conrov et al.l I2006T ). Moreover, at each redshift, 
they have fixed satellite fractions, independent of 
n g . At z = 0, z = 1, and z = 3, our model 
and the observational data cluster near / sat « 0.25, 
/ sat « 0.2, and / sat « 0.1, respectively. We note 
that the lower satellite fractions at high redshift 
are not due to HOD evolution. Figure |4] shows 
that (JVg a t) is higher at high z. Instead, satel- 
lite fractions are lower at high z because all rel- 
evant host halos have A/host > M* and lie on the 
exponentially-decreasing portion of the halo mass 
function, so the relative number of Mi-mass host 
halos to M m i n -mass halos decreases with redshift. 
Nevertheless, these satellite fractions at high red- 
shift are too high to support a power-law galaxy 
correlation function. 

9. Figure [7] implies that the physical mechanisms that 
dictate the HOD of galaxies operate to maintain a 
and / sa t approximately fixed and not to achieve a 
power-law correlation function. 

We have established that the observed power-law correla- 
tion function at low masses and low redshifts should not 
persist at higher masses or redshifts for simple, phys- 
ical reasons. However, exploring the HOD parameter 
space has not revealed a single simple property that 
yields a power-law shape for £(r). In an effort to bet- 
ter understand the factors that drive a power-law £ (r) at 
high precision, we continue to explore the HOD param- 
eter space in a different way. Specifically, we investigate 
the two mass scales in the standard HOD models, M m i n 
and Mi , relative to the characteristic non-linear collapse 
mass, M*. To complement our previous analysis and 
to be consistent with gross theoretical predictions, we 
fix a — 1 and take our two parameters to be M m i n /M* 
and Mi/M m j n . The first ratio specifies roughly the host 
masses that galaxies occupy relative to the exponential 
regime of the halo mass function, and the second ratio 
sets the length of the "plateau" in the HOD. 

Figure probes the power-law £ (r) space as a func- 
tion of the ratios M m i n /M* and Mi/M m - m . For this 
analysis we switch from fixed number density samples to 
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Fig. 8. — Exploration of the HOD parameter space that yields a 
power- law £(r). Each panel corresponds to a different mass thresh- 
old Af m i n (in units of h~ 1 Mq). The y-axis shows redshift (right- 
hand side), which also corresponds directly to the ratio M m ; n /M* 
(left-hand side), since the characteristic non-linear mass M* de- 
pends directly on redshift. The horizontal axis shows the ratio 
Mi/M m i n . We fix the slope of the satellite mean occupation func- 
tion to b e a = 1 and we set the fourth HOD parameter Mq using 
equation 1101 Each point on the horizontal axis therefore corre- 
sponds to a specific set of HOD parameters, while moving along 
the vertical axis shifts the HOD to different redshifts. As in Fig. [7] 
shaded contours represent the 68.3%, 95% and 99.6% power-law 
likelihood spaces and thin solid curves show contours of constant 
satellite fraction. The horizontal and vertical dotted lines corre- 
spond to fixed values of M m i n /M* = 0.05 and Mi/M m - lrL = 30, 
which bisect the best-fit power-law space in all four panels. Solid 
black horizontal lines denote z = 0, below which the parameter 
space corresponds to future epochs. 



fixed mass thresholds, and we choose four values of min- 
imum mass that correspond to a range of sub-L» galax- 
ies (M min = 10 100 , 10 10 - 5 , 10 110 , lO 11 - 5 /^^©), showing 
results for each in a distinct panel. In each panel, we 
sample redshifts from z — —0.9 to z = 2.9 in steps of 
Az = 0.08 (labeled on the right vertical axis). Each red- 
shift value also corresponds to a M m i n /M» ratio, which 
we label on the left vertical axis. At each redshift, we 
also loop over Mi/M m ; n ratios from 1 to 100 in steps 
of A(Mi/M min ) = 2. For every pair of M min /M* and 
Mi/M min values, we compute £(r) using the halo model 
and fit a power-law function in the same fashion as de- 
scribed previously. As before, we show the 68.3%, 95% 
and 99.6% likelihood regions of £(r) consistent with a 
power law (green, blue, and red contours, respectively). 
Also, as before, we show contours of constant satellite 
fraction, / sat (solid black curves). The thick horizontal 
lines at z = are meant to emphasize that the param- 
eter space lying below these lines corresponds to future 
epochs. 

We have repeated this analysis for 
higher mass thresholds (values of M m ; n = 
10 12 , 10 12 - 5 , 10 13 - , lO^ft-^o). However, we do 
not show those results because we found no parameter 
combinations within the 99.6% power-law likelihood 
space. Figure [7] showed that the power-law parameter 
space drifted to higher values of a for lower number 
density (and hence higher mass) samples in an effort 
to drive up the two-halo term to meet the enhanced 
one-halo term. Therefore, it is not surprising that we 
do not find this space when we restrict the slope to be 
a = 1. 

Many interesting results can be drawn from Figure |8l 
Again, we itemize them for the sake of clarity. 

1. In order for £(r) to have a shape consistent with a 
power-law assuming ~ 10% measurement errors, it 
appears necessary for the "plateau" in the HOD 
to be sufficiently long. At all masses and red- 
shifts, Mi/M m ; n needs to be at least ~ 20, oth- 
erwise the large satellite fraction drives a one-halo 
term that is too large relative to the two-halo term. 
Moreover, for past epochs, z > 0, the maximum 
plateau length is Mi/M m - m < 40. Higher values of 
Mi/M m ; n yield a one-halo term that is too weak. 
In fact, Afi/A/ m j n ~ 30 seems to be the preferred 
value to yield a nearly power-law correlation func- 
tion at all masses so long as z > 0. This value 
is denoted by the vertical dotted lines in all the 
panels. 

2. For z > 0, a near power-law £(r) seems to re- 
quire a restricted range of M m i n /M*. Interestingly, 
the value M m i n /M» ~ 0.05 can yield a power-law 
correlation function at all masses for appropriate 
choices of redshift. This value is denoted by the 
horizontal dotted lines in all the panels. This re- 
striction on M m i n /M* means that higher redshift 
samples (when is significantly smaller than to- 
day) can only exhibit power-law behavior if the 
relevant host halos are significantly smaller. This 
possibility becomes irrelevant in a practical sense 
because star formation is inef ficient in small ha- 
los (M < lO 11 /!" 1 ^, e.g.. IConrov fc Wechslen 
[20091: iBehroozi et all l2010YlGuo et all l2010lV )~so 
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they cannot host galaxies that are easily observable 
at high redshift. Figure [8] shows that the lowest- 
mass samples that we consider (top two panels) 
have a nearly power-law £(r) at 1 < z < 2, whereas 
the highest-mass samples have a nearly power-law 
£(r) only at low redshift. 

3. At sufficiently high redshift, near power-law clus- 
tering is no longer achievable at any mass threshold 
corresponding to relatively bright galaxies. Our re- 
sults generally indicate that power-law clustering 
at high redshift can only be achieved if galaxies 
at high redshift occupy halos in a markedly differ- 
ent and more complicated manner than their low-z 
counterparts. 

4. For future epochs these broad results no longer 
hold. A broader range of Mi/M m j n values can be 
made approximately consistent with a power law at 
low values of M min /Af*, or low/negative redshifts. 
For the lowest M m j n samples the power-law likeli- 
hood space is clearly bimodal, with possible ways 
to achieve a power law both at high redshifts and 
at low/future redshifts. 

5. At all masses and redshifts we find that the power- 
law likelihood parameter space has satellite frac- 
tions in the range / sat ~ 0.1 — 0.25, with the 
/sat = 0.15 contour slicing through all of the la re- 
gions. / sa t is naturally strongly dependent on both 
Mi/M m j n and M m i n /M*. Increasing the length 
of the HOD plateau at fixed M m - m and redshift 
makes / sat decrease, as does boosting M m i n /M* 
while keeping M m - m and the plateau fixed. If we 
keep both ratios fixed (i.e., both the HOD shape 
and its position relative to the mass function) then 
the satellite fraction is also approximately fixed, 
regardless of M min . 

7. DISCUSSION AND PRIMARY CONCLUSIONS 

It has been recognized for decades that the two-point 
correlation function has a simple, power-law form with 
£(r) ~ r~ 2 . Observational determinations of galaxy 
two-point clustering spanning more than thirty years 
all yielded results consistent with a single power law 
extending from linear and quasi-linear length scales 
(r > 30/i _1 Mpc) to deeply non-linear scales (r < 
0.1/i _1 Mpc). In this paper, we cast the problem in 
the contemporary setting in which galaxies form in halos 
and subhalos of dark matter and set out to understand 
the physical processes that drive this surprisingly sim- 
ple result. Our primary conclusion is that the nearly 
power-law correlation function of relatively common, L* 
and sub-L* Galaxies at z ~ is a coincidence and does 
not reflect any general principle of structure formation 
or galaxy evolution. So how did we arrive at this conclu- 
sion? 

First, the efficiency of galaxy formation is dependent 
upon halo mass and it has been determined both theoret- 
ically and empirically that there is a halo mass scale be- 
low which galaxy formation is inefficient, roughly Af ga i ~ 
1Q io.5 h -i Mp) (IConroy fc Wechslerll200l IBehroozi et all 
[2(jT0tlGuo et^l20lQlh A number of things can set this 
scale including atomic and molecular physics and feed- 
back from supernovae and active galactic nuclei (for a 



recent review article see lBenso"nll2010D . This mass scale 
is Mgai < M*, so L* and sub-L* galaxies are common. 
Had Mgai been greater than or similar to M* , most bright 
galaxies would lie in comparably rare halos and be rare 
themselves. In such a case, one-halo clustering would be 
too strong to be compatible with a power law. M* is 
not determined by galaxy formation physics but is set by 
the completely unrelated processes that establish the am- 
plitude of cosmological density fluctuations, presumably 
primordial inflation. 

Second, power-law clustering requires that some of the 
galaxies formed within relatively large subhalos are de- 
stroyed. Destruction is due primarily to mass loss, and, 
to a lesser extent, merging with the central galaxy as 
a result of dynamical friction. Without this destruction, 
satellite fractions would be too high and small-scale clus- 
tering too strong compared with large-scale clustering. 
In a forthcoming paper, we perform more sophisticated 
modeling to make the connection between subhalo mass 
loss and stellar mass loss in order to make predictions for 
the amount of intracluster light. Large-scale clustering 
is principally set by large-scale matter density fluctua- 
tions and is insensitive to the details of galaxy formation 
within halos, while the strength of small scale clustering 
grows in proportion to the fraction of galaxies that are 
satellites and in inverse proportion to the number den- 
sity of the galaxies of interest. As it turns out, precisely 
the right amount of subhalo destruction has occurred by 
redshift z ~ in a concordance cosmology to produce a 
single, unbroken, power-law £(r). 

Evolution of the satellite fraction is set by a compe- 
tition between halo mergers, which increase / sa t, and 
destruction by dynamical processes, which occur on a 
dynamical timescale and reduce / sat . At high redshifts, 
mergers occur more rapidly than destruction for halos 
with masses > Af ga i. The low-redshift merger rate de- 
clines in part due to the fact that M ga i < Af* at z < 1. 
Halos with masses below M* become relatively more 
likely to merge with a larger object than to acquire 
new substructure co mpared to coun terparts with masses 
greater than M* (see IZentnerl 120071 ). More importantly, 
the rate of halo mergers is quenched at z < 1 as dark 
energy begins to suppress further cosmological structure 
growth. As merger rates decline, satellites are depleted 
with time. Therefore, at z ~ 0, the correlation function 
is nearly a power law because the competition between 
the accretion and destruction rates has struck just the 
right balance to yield the appropriate value of / sa t . 

The merger and destruction rates will once again be- 
come unbalanced in the future as halo merging is stifled 
by dark energy and existing satellite galaxies are slowly 
destroyed over many dynamical times through complex 
interactions in their host environments. We show that 
this will result in small-scale clustering that will be sig- 
nificantly too weak to be consistent with a power law. 

Largely as a consequence of the merger/destruction 
competition, £(r) evolves through cosmic time, achiev- 
ing a power law only near z ~ for L ~ L* and dimmer 
galaxies. The processes of galaxy formation, the am- 
plitude of cosmological density fluctuations, the abun- 
dance of dark matter, and the nature of the dark energy 
are thought to be completely distinct and determined by 
unrelated physics. So the power-law £(r) at z ~ is a 
coincidental conspiracy. 
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In establishing these broad conclusions, we have per- 
formed an exhaustive investigation of the ingredients of 
the galaxy correlation function, which has revealed many 
interesting, more detailed conclusions. These can be 
summarized as follows. 

1. We find that satellite halo mass loss is the princi- 
ple dynamical process responsible for depleting suf- 
ficient substructure so as to nearly align the one- 
and two-halo terms to yield a power-law correla- 
tion function at low redshift. Dynamical friction 
plays a smaller supporting role, accounting for an 
additional ~ 15% of subhalo destruction. 

2. The shape of the correlation function is strongly 
mass dependent. For instance, at low redshift de- 
viations from a power law £(r) grow with increas- 
ing host halo mass. This drives stronger devia- 
tions from a power law for higher luminosity galaxy 
samples. The best power-law fits derived from our 
model are for galaxies residing in halos that are 
common enough to correspond to ~ L* and dim- 
mer galaxies, in agreement with observations. 

3. The correlation function is highly redshift- 
dcpcndcnt. The sensitivity of the one- halo term 
to the HOD, coupled with the relative insensitiv- 
ity of the two-halo term, implies that achieving a 
power-law requires fine-tuning the number of satel- 
lite galaxies per halo. The satellite galaxy abun- 
dance evolves with redshift, driven by the evolving 
balance between accretion and destruction, with an 
enhanced amount of substructure at high redshift. 
Therefore, the correlation function can only achieve 
a power law during those epochs when substructure 
has evolved to align the one- and two-halo terms. 
The correlation function is boosted on small scales 
at high z, the one- and two-halo terms join at z = 
to form a power-law, then the power law is once 
again broken in future epochs. 

4. For three chosen number densities corresponding to 
low-redshift, ~ L* and dimmer galaxies, we probed 
the most likely power-law space as a function of 
redshift for a parametrized HOD. We find that 
there is a relatively narrow range of satellite frac- 
tions for £(r) to be consistent with a single power 
law (assuming ~ 10% measurement errors) at any 
given redshift. At all redshifts and masses, power- 
law correlation functions have satellite fractions in 
the range / sat <~ 0.1 — 0.25. It is difficult to achieve 
a power-law correlation function at z > 3 for any 
number density. 

5. We find that to achieve a power law £(r) at high 
mass or redshift, the slope a of the satellite galaxy 
occupation function must be significantly steeper 
than unity (for instance, greater than 2 at z = 3). 
This would imply that the mapping of galaxies to 
halos is much more complicated than we think, 



since the number of galaxies would have to be very 
different than the number of subhalos of a particu- 
lar size. Instead, it appears that the processes that 
govern galaxy formation do not care about the con- 
ditions needed to achieve a power law £(r). 
6. The ratio Mi/M min (the "plateau" of the HOD) is 
a key ingredient for predicting the shape of £(r). 
The prominence of the plateau is a measure of sub- 
structure abundance. Along with Mi/M m i n , it is 
also necessary to characterize the ratio M m i n /M*, 
which specifies what halo masses galaxies occupy 
relative to the halo mass function. By main- 
taining the combination of Mi/M m - m <~ 30 and 
A^min/^* ~ 0.05 we can achieve a near power law 
for redshifts in the range 0—1.5 and the appro- 
priate mass threshold at each redshift (the mass 
threshold is M min ~ M*/20, with M* set by the 
redshift). At higher redshifts this criterion is met 
for galaxies that are most likely too dim to be 
observed. For example, achieving the requisite 
M m j n ~ M*/20 at z — 2 corresponds to a halo 
mass of M m j n ~ 10 9 /i _1 M in which star forma- 
tion is inefficient. 

This work has allowed us to formulate a general picture 
of the nature of the galaxy two-point correlation func- 
tion. Halo abundances and subhalo populations evolve 
with time. At high redshifts, halos large enough to har- 
bor galaxies are rare and subhalos are abundant within 
these hosts. With time, host halos that harbor galaxies 
generally become more common (though the specifics of 
this evolution can be subtle) and subhalos within these 
hosts become relatively less abundant. All the while, 
large-scale matter correlations grow, but the clustering 
bias of large halos evolves to largely compensate for this 
large-scale growth of structure. These effects, considered 
either individually or in tandem, change the HOD and 
the shape of £(r). As a result, the correlation function 
evolves through an epoch where it is close to a power 
law and this epoch happens to be near z ~ 0. From 
our broad discussion and detailed conclusions, it is clear 
that a nearly power-law correlation function requires a 
conspiracy between otherwise unrelated processes such 
as the early Universe physics that established the initial 
conditions for low redshift structure, the detailed phys- 
ical processes that determine galaxy and star formation 
efficiency, and the growth rate of cosmic structure set 
largely by the abundances of dark matter and dark en- 
ergy. The low-redshift power-law galaxy two-point func- 
tion is thus a mere cosmic coincidence. 

We appreciate many helpful discussions from James 
Bullock, Scott Drake, Andrcy Kravtsov, Cameron 
McBride, Jeff Newman, Risa Wechsler, and Zheng 
Zheng. AAB is supported by Vanderbilt University and 
the Alfred P. Sloan Foundation. The work of ARZ is 
funded by the University of Pittsburgh and by the Na- 
tional Science Foundation through grant AST 0806367. 
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